clear;

Hd = filter1;

f0=100e6;
T = 1 / (10000*f0);
t = 0:T:(15 * 1/f0);
lt = length(t);


stdn=sqrt(13);
N=stdn*randn(1,lt);

y = filter(Hd, N);

fft_sh = fft(N);
X1 = abs(fft_sh);
X = X1/max(X1);

fft_s = fft(y);
Y1 = abs(fft_s);
Y = Y1/max(Y1);
%Sn = fft(N);
%S = fft(y);

f = 0:f0/15:1/T;
%plot([Sn; S]')

plot (f,20*log10([X; Y])')